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TRANSONIC  FLOW 

FLOWFI  ELD  OSCILLATIONS 

TWO-DIMENSIONAL  MODEL 

AIR  BREATHING  PROPULSION  SYSTEMS 

INLET  FLOW  DISTORTION 
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The  second  «nd  third  veers  of  the  contract  comprise  Its  second  phase,  aimed  at 
exploring  the  effects  of  periodic,  downstream  excitation  on  the  transonic  flow 
a  two-dimensional  diffuser  model.  The  present  report  cover*  the  second  year  of 
the  contract.  The  pulse  generator  was  incorporated  In  the  diffuser  model,  and 
Its  controls  were  synchronised  with  the  optical  lnat rumentat Ion.  The  boundary- 
layer  control  lystm  was  turned  for  best  two-dimensionality  of  the  flow.  The 
flowfleld  Instrumentation  was  selected,  and  an  appropriate  actuator  was 
constructed.  Perturbation  of  the  flow  by  the  actuator  was  determined  to  be  - 
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acceptably  small.  Boundary- layer  profiles  were  determined  on  all  four  walla, 
and  the  Hach  number  distribution  was  mapped  In  detail  over  the  diffuser  exit 
cross-section.  Surface  pressure  distributions  were  measured,  and  spark 
schlleren  photos  were  taken  for  the  available  pressure  ratio  range. p. 
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The  work  reported  herein  was  conducted  by  the  McDonnell  Douglaa  Research 
Laboratories  (MDRL) ,  St.  Louis,  Missouri,  for  the  United  States  Air  Force 
Office  of  Scientific  Research,  Bolling  Air  Force  Base,  D.C.,  under  contract 
F49620-77-C-0082.  The  *»rk  reported  was  conducted  from  1  April  1978  to 
31  March  1979  In  the  Flight  Sciences  Department,  managed  by  Dr.  R.  J.  Hakklnen. 
The  principal  Investigator  was  Dr.  M.  Sa  l  bon,  and  the  co- Invest  lgator  was 
Mr.  J.  C.  Kroutll.  The  program  technical  manager  was  Dr.  D.  Samaras,  Air  Force 
Office  of  Scientific  Research. 

This  technical  report  has  been  reviewed  and  Is  approved. 


£  j?. 

R.  J.  hakklnen 

Chief  Scientist,  Flight  Sciences 
McDonnell  Douglas  Research  laboratories 


D.  P.  Ames 

Staff  Vice  President 

McDonnell  Douglas  Research  Laboratories 


' 


yCr 


'L*>C  Tib 
Vi  ■  . 


■u 


unc-d 

ls-Ul  oa 


By 


~  -  _ _ 

r<’Cc.ai 


fl 


111 


table  of  contents 

P>M 

1.  INTRODUCTION .  1 

2.  RESEARCH  OBJECTIVES  .  2 

3.  RESEARCH  STATUS  .  3 

3.1  Diffuser  Model .  3 

3.2  Perturbation  Generator  .  •  .  3 

3.3  Flovfield  Diagnostic  Method .  *> 

3.4  Boundary  Layer  Control  (BLC)  Tuning  .  2 

3.5  Baseline  Data  -  No  Excitation .  2 

3.6  Summary .  21 

4.  PUBLICATIONS .  23 

5.  PROFESSIONAL  PERSONNEL  .  14 

6.  INTERACTIONS .  23 

7.  NEW  DISCOVERIES .  26 

8.  RELATED  RESEARCH  .  17 

9.  REFERENCES .  18 

10.  DISTRIBUTION  LIST .  19 


V 


LIST  OF  ILLUSTRATIONS 


Figure  Page 

1  Dimensions  of  the  diffuser  model  .  3 

2  Diffuser  model  Installed  In  the  Internal  Fluid  Dynamics 

Facility  .  4 

3  Disturbance  generator  used  to  periodically  modulate 

diffuser  exit  area .  S 

4  Top-wall  surface-pressure  distributions  for  various 

pressure  ratios  .  8 

5  Total-pressure  profile  of  boundary  layer  at  top  inlet 

center  .  9 

t>  Total-pressure  contour  map  for  exit  station  at  V  •  1.205  .  .  9 

7  Mean  shock  position  and  shock  amplitudes  as  functions  of 

pressure  ratio  .  10 

8  Wall-pressure  fluctuation  Intensity  (rms)  at  various  stremwlse 

locations,  (u  *  local  mean  velocity) .  11 


vl 


1.  INTRODUCTION 


The  present  study  Is  part  of  a  continuing  Investigation  of  unsteady, 

transonic  diffuser  flows,  with  application  to  dynamic  distortion  In  the 

Inlets  of  fighter  aircraft  and  airbreathing  missile  propulsion  systems. 

The  Investigation  Is  focused  on  a  simple,  two-dimensional  diffuser 

configuration  displaying  a  weak  shock  shortly  downstream  of  the  throat. 

Flows  In  this  configuration  were  found  to  display  self -excited  oscillations 

Involving  the  shock  and  the  entire  subsonic  flow  behind  It.  The  shock 

displacement  amplitudes  are  comparable  to  the  throat  height  and  occur  at 

frequencies  character  lied  by  the  flow  speed  and  the  length  of  the  divergent 

diffuser  section.  The  time-mean  and  fluctuating  properties  of  this  flow  and 

1-4 

the  dependence  of  these  properties  on  shock  strength  have  been  documented 

The  first  year  of  the  contract  was  devoted  to  exploring  the  effects  of 
approach  boundary-layer  thickness.  This  work  Is  complete  and  has  been 
partially  documented. 

The  second  and  third  years  constitute  the  second  ptiase  of  the  contract. 
The  goal  of  this  phase  Is  to  explore  the  response  of  the  diffuser  to  t  lmc- 
dependent,  downstream  boundary  conditions.  The  primary  purpose  Is  to  provide 
benclwark  data  for  testing  related  theoretical  efforts  In  progress  at  MDRI.  and 
elsewhere.  The  experimental  results  should  also  help  define  the  role  of  the 
Inlet  diffuser  response  In  coupled  lnl  ct/combu  stor  oscillations  In  ramjets*1. 


2.  RESEARCH  OBJECTIVES 


The  specific  objectives  of  the  second  phase  of  the  contract  are  as 
follows: 

1)  Design  and  build  appropriate  hardware  suitable  for  creating 
nearly  planar,  periodic  waves  downstream  of  the  diffuser  model. 

2)  Obtain  detailed  steady  and  dynamic  measurement s  of  the  surface 
pressures  and  shock  positions  at  various  pressure  ratios  and 
perturbation  parameters. 

3)  Obtain  phase-averaged  flowfleld  velocity  and  velocity  fluctuation 
measurements  at  selected  test  parameter  coeiblnat Ions. 

4)  Reduce,  analyze,  and  document  data  In  a  form  useful  for  comparison 
with  theoretical  prediction  methods. 


3.  RESEARCH  STATUS 


3_.  1  J)lffuiier  Mod  e  1 

The  diffuser  model  used  In  this  research  (Figure  1)  was  available  from 
a  parallel,  MDRL-sponsored  experimental  effort.  The  model,  being  a  third- 
generation  design.  Incorporates  the  experience  from  three  years  of  similar 
work.  The  ex it-to-t hr oat  area  ratio  is  1.52,  the  ratio  of  divergent  length 
to  throat  height  Is  7.2,  and  the  aspect  ratio  of  the  throat  cross-section 
Is  4 .  Throat  height  Is  44  m».  Sidewall  and  bottom  boundary-layer  growth  is 
limited  by  three  sets  ol  ram-type  suction  slots,  with  Independently  controlled 
and  measured  flow  rates  through  each  set.  The  sidewalls  are  2. 54-cm-thlck, 
sc  hi  ieren-qun 1  It  y  glass  panels. 
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Figure  t.  Dimeruion*  of  the  cliff  u»er  model. 
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The  model  Is  supplied  with  clean,  dry  air  through  a  plenum  cluunber 
equipped  with  a  25:1  area  ratio  contraction  section.  The  air  exhausts  to  the 
laboratory,  keeping  the  exit  pressure  constant.  Shock  strength  is  varied  by 
controlling  the  plenum  chamber  pressure  and  thereby  the  overall  pressure 
ratio,  V  (ratio  of  plenum  total  to  exit  static  pressures).  Shock  strength 
increases  raonotonlcally  with  v;  the  strongest  obtainable  shock  Mach  number  is 

msh  "  l-35‘ 

The  model  is  extensively  instrumented  with  numerous  surface-pressure 
measuring  orifices  and  hlgh-rcsponse  surface-pressure  transducers.  An  optical 
system  with  two  51-cm  d  lam  concave  mirrors  Is  used  to  obtain  shadovgjaph  and 
schlleren  photographs,  and  motion  pictures.  Used  in  connection  with  a  line- 
scan  imaging  earner a  and  appropriate  electronics,  the  system  also  provides  a 
time-dependent  analog  signal  proportional  to  streamwlsc  shock  position. 

A  photograph  of  the  model.  Installed  and  fully  Instrumented,  Is  shown 
In  Figure  2. 
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Figure  2.  Diffuser  model  installed  in  the  Internal  Fluid  Dynamics  Facility. 

Pulse  generator  motor  and  encoder  visible  on  near  and  far  sides  of  model,  respectively. 
Ducts  below  model  connect  BLC  suction  slot  manifolds  to  vacuum  pump. 
Components  of  schlieren  system  are  visible. 
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3.2  Perturbation  Generator 


A  detailed  study  of  possible  methods  to  generate  perturbations  led  to 
selection  of  the  mechanical  device  shown  in  Figure  3.  The  active  element  is  a 
rotor  of  rectangular  cross  section,  embedded  in  the  bottom  wall  of  the  model 
and  rotated  around  its  axis  by  a  variable-speed  motor.  The  rotor  modulates 
the  exit  area  from  0  to  -9.5Z,  which  was  found  to  produce  a  4-cm  shock 
d lscplacement  amplitude.  Frequencies  are  continuously  variable  to  150  Hr,  and 
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Ftfura  3.  Disturbance  generator  used  to  periodically  modulate  diffuser  exit 


locked  to  present  a  flat  boundary  surface  to  the  wall,  corresponding  to  the 
baseline  case  of  perturhat Ion  f rec  flow. 

The  area  modulation  produces  pressure  waves  tliat  propagate  In  the 
upstream  direction.  As  the  wave  fronts  recede  from  their  point  of  origin, 
their  radius  of  curvature  Increases  and  they  become  Increasingly  more  planar. 
The  wall  pressures  are  monitored  on  both  top  and  bottom  walls  at  the  exit 
station,  located  approximately  three  duct  heights  upstream  of  the  generator. 

The  pressure  wave  forms  are  recorded,  phase  averaged,  and  used  as  exit 
boundary  conditions  for  related  theoretical  efforts.  Bv  examining  the  phase 
relationship  between  the  lop  and  bottom  probe  signals,  the  planarity  of  the 
wave  front  can  be  established. 

An  optical  encoder  is  mounted  on  the  rotor  shaft  with  dual  outputs:  1 
pulse/revolution  and  360  pul se/revolut Ion  signals  are  available.  The  1  pulse/ 
revolution  output  of  the  encoder  Is  used  as  a  synchronous  signal  to  initiate 
the  sampling  period  In  the  data  reduction  process.  The  360  pu 1 ses/revolut Ion 
output  Is  used  to  define  phase  angles  In  data  sampling.  In  conjunction  with  an 
adjustable  counter  circuit.  The  latter  signal  Is  also  used  to  drive  a  stepping 
motor  mounted  elght-bladed  r!»opper  wheel  placed  on  the  axis  and  near  the  slit 
of  the  schllcren  system.  The  light  beam  Is  periodically  Interrupted  and  the 
schlleren  image  Is  strobed  to  provide  real-time,  slew-motion  observation  of 
the  shock  motion.  This  system  references  strobing  to  the  excitation  frequency 
and  eliminates  problems  caused  by  minor  drifts  in  either  frequency.  The  drive 
motor  of  the  generator  and  part  of  the  rotor  are  shown  in  figure  2,  near  the 
d  1  f fuser  exit. 

3_._3  How-field  Dlagnost  lc  Met  hod 

The  primary  goal  of  this  project  Is  to  explore  the  dependence  of  the  flow 
on  external  parameters  (shock  strength,  excitation  frequency  and  amplitude). 

In  order  to  establish  the  respective  trends,  numerous  test*  must  he  made 
Involving  extensive  variations  of  flow  parameters.  The  demands  of  reliability 
and  speed,  implicit  In  such  a  program,  led  to  the  selection  of  a  miniature 
total-pressure  probe  to  determine  the  core  flow  velocity  distributions 

(u  (x,t))  for  a  number  of  parameter  combinations. 

c 

The  sensor  used  Is  a  Kullte  model  XCQ-062  pressure  transducer,  mounted  on 
a  sting  at  a  fixed  vertical  distance  (19mm)  from  the  bottom  wall.  The 
actuator  Is  huilt  Into  a  removable  panel  fitted  to  the  bottom  wall.  The 


sensor  diameter  la  1.6  tan  and  causes  no  measurable  disturbances  In  the  wall 
pressure  distributions,  with  the  exception  of  an  approximately  2-cm  Ion# 
region  on  the  bottom  wall  behind  the  sting. 

The  probe  signal  Is  recorded  on  I'M  tape  simultaneously  with  reference 
pulses  from  the  encoder  of  the  perturbation  generating  system,  thereby 
establishing  a  common  time  base.  Phase  averages  are  determined  by  processing 
the  d  igit  Iced  data. 


3 . 4  Boundary  Layer  Control  (Bl.C)  Tuning 

In  a  series  of  preliminary  test  runs,  the  flow  rates  removed  through  each 
of  the  three  slot  sets  were  measured  as  functions  of  the  respective  suction 
manifold  pressures  (normalized  by  the  total  pressure  at  the  diffuser  Inlet). 

In  normal  operat  Ion,  the  manifold  pressures  are  monitored  and  the  Individual 
bleed  flow  rates  arc  computed  using  the  calibrated  relations.  The  sum  of 
the  computed  bleed  rates  Is  compared  with  the  total  bleed  flow  which  Is  also 
Independently  monitored.  The  computed  and  measured  flow  rates  routinely  agree 
within  about  S-6T,  which  Is  satisfactory.  This  arrangement  requires  only  one 
flow  meter  for  the  BLC  system  Instead  of  three. 

The  fraction  of  the  flow  removed  at  each  slot  Is  approximately  equal  to 
the  ratio  of  the  slot  area  to  the  cross-sectional  area  of  the  diffuser  model 
at  the  slot.  The  respective  area  fractions  for  the  three  slot  sets  are: 

0.05®,  0.023  and  0. 12B.  No  fluid  Is  removed  between  the  throat  and  the  exit 
station;  thus,  mass  flow  Is  conserved  In  the  experiment,  as  also  assumed  In 
the  related  theories. 

3  .  S  Basel  lne  Pat  a  -  No  F.xc  ltat  ton 

All  data  In  the  second  contract  vear  were  taken  without  excitation  to 
establish  a  baseline  for  the  assestosent  of  response  character  1st  lea. 

A  family  of  normalized  surface  pressure  distributions  Is  shown  In 
Figure  A.  These  data  are  repeatable  with  high  precision  (♦0.002).  The  "bump" 
near  the  throat  Is  caused  by  a  machining  Imperfection  In  the  top  wall  contour. 
The  Irregularity  Is  noticeable  but  Inconsequential,  and  no  attempt  was  made 
to  eliminate  it.  The  increasingly  sharp  "kink"  at  the  beginning  of  the  post¬ 
shock  subsonic  region  Indicates  the  appearance  of  shock- Induced  separation. 

The  experiment  thus  Is  capable  of  producing  both  attached  and  separated 
transonic  flows. 
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Inlet  velocity  profile*  were  measured  using  a  miniature  total  pressure 
probe  (0.08  am  *  0.32  am  frontal  cross  section).  The  profiles  (Figure  5)  are 
well  behaved  and  thin:  the  blockage  at  the  Inlet  station  (*  •  -5.1)  Is 
approximately  2.22. 

The  exit  station  velocity  profiles  were  measured  with  a  total  pressure 
probe  rake  having  11  total  pressure  lubes  (0.81  am  o.d.)  aligned  vertically. 
The  rake  wns  traversed  horizontally  to  cover  the  entire  cross  section 
(Figure  6) .  Tbe  Influence  of  the  sidewall  boundary  layers  on  the  central 
portion  of  the  flow  Is  clearly  minimal.  The  residual  deviation  from  two- 
dimensional  lty  on  the  right  side  of  Figure  6  was  traced  to  an  asymmetry  In  the 
second  suction  slot  and  has  been  corrected. 

Preliminary  shock  position  and  oscillation  amplitude  data  were  obtained 
by  visual  observation  of  the  schlieren  image  of  the  flow  (Figure  7).  For 
weak  shocks,  the  flow  Is  remarkably  stable,  and  the  shock  displays  minor 
positional  deviations  only.  As  the  shock  strength  Is  Increased,  Increasingly 
larger  amplitude  oscillations  develop,  but  even  the  largest  amplitudes  (at 
V  "  1.36)  are  only  approx  lmately  ♦  101  of  the  throat  height.  This  amplitude 
Is  much  smaller  than  that  Introduced  by  the  perturbation  generator,  which 
should  facilitate  separating  of  the  two  types  of  motion  in  the  data  analysis. 
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Figur*  6.  Total  pranura  contour  map  for  axit  station  at  »•  -  1 .206. 
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Figure  7.  Mean  shock  position  and  shock  amplitudes  as  (unctions  of  pressure  ratio. 

Selected  fluctuating  surface  pressure  (nns)  Intensities  are  shown  In 
Figure  8.  The  Inlet  pressure  fluctuations  are  constant  at  approximately  IX 
of  the  local  dynamic  pressure.  Fluctuations  at  the  exit  stations  Increase 
with  pressure  ratio;  steep  growth  begins  at  v  *  1.22  (which  also  marks  the 
appearance  of  appreciable  shock  oscillations). 

A  sensor  (labeled  exhaust  in  Figure  8)  was  placed  25.4  on  upstream  of  the 
channel  end.  The  fluctuation  level  at  this  location  Is  high.  Indicating  that 
the  conventional  assumption  of  equating  exhaust-station  static  preamire  to 
ambient  pressure  may  he  significantly  violated  In  this  experiment. 


Figure  8.  Wall  prmtura  fluctuation  intensity  (rim)  at  various  streamwise  location*. 
(0  •  local  mean  velocity). 


3 . 6  Summary 

All  necessary  equipment  and  lnatrumentat ion  required  to  achieve  the  con 
tract  objective*  ha*  been  built,  assembled,  checked  out,  and  placed  In 
operat Ion. 

A  large  portion  of  the  necessary  baseline  Information  has  been  acquired 
Data  acquisition  Is  In  progress. 


The  exper  Imental  configuration  was  demonstrated  to  possess  two  features 
highly  favorable  to  achieving  the  program  objectives.  First,  the  natural 
fluctuations  have  low  amplitudes  and  thus  will  not  Interfere  with  the  study  of 
externally  Introduced  oscillations.  Second,  flows  both  with  and  without 
separation  can  be  produced,  allowing  the  assessment  of  effects  attributable 
to  flow  separation. 
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4.  PLICATIONS 


A  sinuitcr  lpt  describing  the  result*  of  the  first  phase  of  the  contract 
Is  in  preparation.  Thla  paper  la  Intended  for  the  A1AA  Journal  and  la 
t  it  led : 

a  The  Effects  of  Approach  Boundary  Layer  Propertlea  on  Shock 

Oac illations  in  a  Diffuser  Flow,  by  H.  Sajben  and  J.  C.  Kroutll. 
Related  research  work,  sponsored  by  MDRL,  led  to  two  papers,  both 
accepted  by  the  A1AA  Journal: 

•  Real-Time  Optical  Measurements  of  Time-Dependent  Shock  Position, 
by  M.  Sajben  and  R.  C.  Crltea 

a  Shock-Wave  Oscillations  in  a  Transonic  Diffuser  Flow,  by  C.  P.  Chen, 
M.  Sajben,  and  J.  C.  Kroutll. 
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6.  INTERACTIONS 


Several  seminars  were  given  by  M.  Sajben,  which  Included  material 
resulting  from  the  present  contract. 

1)  USAF,  Flight  Dynamics  Laboratory 

Wr ight-Patterson  Air  Force  Base,  Dayton,  OH 

Host:  Dr.  G.  K.  Richey 

19  March  1979 

2)  University  of  Michigan 

Ann  Arbor,  MI 

Host:  Prof.  T.  C.  Adamson 

20  March  1979 

3)  NASA-Lewis  Research  Center 

Cleveland,  OH 

Host:  Dr.  B.  H.  Anderson 

22  March  1979 

4)  Case  Western  Reserve  University 

Cleveland,  OH 

Host:  Professor  I.  Grcber 

23  March  1979 

Experimental  results  f rom  the  Initial  stages  of  this  work  have  been  used 
to  validate  numerical  ccimputat Ions  carried  out  at  NASA-Ames  Research  Center. ^ 
The  work  at  Ames  was  done  by  Dr.  T.  J.  Coakloy,  who  used  several  types  of 
turbulence  models  In  time-averaged  Nav ler-Stokes  calculations  to  predict  flows 
for  the  MDRL  diffuser  model. ^  The  predicted  time-mean  flows  agree  well  with 
the  experimental  data. 

The  experience  and  understanding  gained  from  this  work  enabled  M.  Sajben 
to  provide  useful  Inputs  to  McDonnell  Douglas  Astronautics  Company  (MDAC)  In 
connection  with  the  design  of  a  ram-jet  missile  propulsion  system.  The 
ram-jets  considered  by  MDAC  employ  dump-combustors  which  are  prone  to  develop 
screech  or  burr.**  It  Is  believed  that  the  coupling  between  the  combustor  and 
the  diffuser  is  an  Important  factor  In  determining  this  type  of  Instability. 
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7.  SEW  DISCOVERIES 


So  Invention  disclosures  were  filed  during  the  second  year  of  this 
contract.  The  results  obtained  are  of  a  scientific  import  with  a  primary 
purpose  to  assist  development  of  predictive  theories. 
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8.  RELATED  RESEARCH 


The  contract  work  la  complemented  by  a  parallel  IRAD  effort,  which  was 

Initiated  In  1976.  During  the  second  year  of  the  present  contract,  this  IRAD 

effort  Included  both  experimental  and  theoretical  contributions. 

The  experimental  work  dealt  with  the  effects  of  area  ratio  on  the 

dynamics  of  diffuser  flow  and  established  the  generally  destabilizing 

effects  of  increasing  area  ratio.  The  Improvement  of  velocity-field  measure- 

8 

sent  capabilities  at  MDRL  also  was  an  objective. 

The  theoretical  effort  (by  Dr.  Meng-Slng  Llou)  produced  an  operational 
computer  code  suitable  for  predicting  time-dependent,  viscous,  transonic  flow 
In  a  two-dimensional  dlffusgr,  subject  to  time-dependent  boundary  conditions. 

g 

The  theory  Is  capable  of  predicting  buzz-like  instabilities. 
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